Projection of plane-wave calculations into atomic orbitals 
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The projection of the eigenfunctions obtained in standard plane-wave first-principle electronic- 
structure calculations into atomic-orbital basis sets is proposed as a formal and practical link between 
the methods based on plane waves and the ones based on atomic orbitals. Given a candidate atomic 
basis, (i) its quality is evaluated by its projection into the plane-wave eigenfunctions, (li) it is 
optimized by maximizing that projection, [lii) the associated tight-binding Hamiltonian and energy 
bands are obtained, and [iv) population analysis is performed in a natural way. The proposed 
method replaces the traditional trial-and-error procedures of finding appropriate atomic bases and 
the fitting of bands to obtain tight-binding Hamiltonians. Test calculations of some zincblende 
semiconductors are presented. 
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Plane waves and atomic orbitals provide two widely 
used methods for electronic structure calculations, each 
of them showing advantages and disadvantages. This pa- 
per proposes a way to use the advantages of each method 
to complement the other. 

On one hand, plane waves (PW), in addition to ah mi-. 
tio pseudopotentialsEl provide a very successful schemecl 
to calculate the ground state properties of a large num- 
ber of systems giving accurate results for electron densi- 
ties, total energies, atomic forces, geometries, and energy 
bands. The origin of its wide applicability is the flexibil- 
ity, simplicity, and accuracy of the plane- wave basis set: 
it has a homogeneous and universal completeness con- 
trolled by a single parameter, the energy cut-off. It is a 
floating basis, independent of atomic positions, and this 
greatly facilitates the calculation of atomic forces. Being 
equally complete everywhere, however, makes the PW 
basis inefficient, requiring a relatively large basis size for a 
given accuracy. In addition, plane waves require periodic 
boundary conditions. When the system is not periodic 
(molecules, surfaces, defects) such a periodicity must be 
imposed artificially by the use of supercells, what further 
increases the inefficiency of the basis. 

On the other hand, methods based on the linear com- 
bination of atomic orbitals (LCAO) are far more effi- 
cient in terms of basis size, since atomic orbitals are 
much better siuted to represent the molecular or Bloch 
wavefunctions.u In addition, LCAO bases have recently 
experienced a renewed interest because of their suitabil- 
ity for order-N methodsa in which computational effort 
scales linearly with system size. Local bases also of- 
fer a natural way of quantifying magnitudes like atomic 
charge, orbital population, bond charge, charge transfer, 
etc., through population analysis. El These quantities are 
essential for the chemical analysis of a condensed matter 
system since most of the chemical language is based on 
them. They provide a deep physical insight. 

The price to be paid for the advantages mentioned 
above is the greater difficulty in the choice of the LCAO 
basis: since the results of calculations are affected by the 



election of the basis, a testing of the completeness and 
quality of the basis is critical. This task is much worse 
defined than in the case of plane waves since there are 
many more parameters than just a plane- wave energy 
cutoff, and usually ends up in a very cumbersome and 
non systematic trial-and-error procedure. 

Thus, due to the difficulties and ambiguities of the 
LCAO basis construction and their less well controlled 
completeness, they are frequently considered more ap- 
proximate than plane waves, and they have been tradi- 
tionally used in less accurate tight-binding type of calcu- 
lations. However, once a well optimized LCAO basis has 
been obtained, there is no reason why highly accurate 
ah initio calculations cannot be performed with it.Q This 
rises the question of whether plane- wave calculations can 
be used to find an optimal or near-optimal LCAO basis, 
which can then be used in more complex situations while 
keeping a feasible computational effort. 

Atomic orbitals constitute non-orthogonal basis sets in 
molecules and solids. Even though this kind of bases are 
more cumbersome to handle than orthogonal ones, they 
have shown to be more transferable and less environment- 
dependent, and to have a shorter range of interactions, 
being thu&_more adequate for non-parametrized LCAO 
methods.Lfp Non-orthogonalized atomic basis sets are 
therefore used throughout this work.cl 

Most schemes proposed to optimize LCAO local bases 
are based on two procedures: (i) the minimization of to- 
tal energies in atoms, molecules, or solids,E£l and {it} the 
optimization of energy bands as compased with exper- 
iments or with plane-wave calculations.Enlj In contrast, 
the present work links the plane-wave and LCAO meth- 
ods by projecting the eigenstates of a plane-wave cal- 
culation into the Hilbert space spanned by the atomic 
basis. This provides an efficient and practical analyti- 
cal tool which is applied to (?) characterize and optimize 
atomic basis sets, (««) obtain LCAO bands and Hamilto- 
nians, and iiii) analyse electronic structures by means of 
population analysis.113 

In this work the reference plane-wave calculations are 
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performed within the local density appcpximation (LDA) 
for electron exchange and correlaticpjES using pseudopo- 
tentials that replace core electronsJlil These approxima- 
tions are, however, not essential to the method. They can 
be replaced by others at will, the only key factor being 
the expansion of the one-particle wave functions in plane 
waves. The theory presented in this paper considers the 
PW reference calculations to be converged in plane-wave 
energy cutoff. The explicit consideration of the uncom- 
pleteness of the PW basis will be the subject of a future 
study. 

Given the PW results for the system of interest, and 
given the set of Hamiltonian eigenstates to be considered 
for the projection, the quality of an atomic basis is quan- 
tified by its ability to represent those eigenstates, i.e., 
by measuring how much of the subspace of the Hamilto- 
nian eigenstates falls outside the subspace spanned by the 
atomic basis. For that we define the spillingiB parameter 
S 

^ = ]^ aT E E(^"WI(1 - ^(k))IV'o(k)) (1) 

k a—1 

where [^/'^(k)) are the PW calculated eigenstates, and 
Nk and Na are the number of calculated k points in 
the Brillouin zone and the number of bands considered, 
respectively. P(k.) is the projector operator into the sub- 
space of Bloch functions of wave vector k, generated by 
the atomic basis, and defined as usual for non-orthogonal 
basis 

p(k) = J2 - E \M^))s-.'ik){M^)\, 

(2) 

where 

V(k) = (0^(k)|0,(k)) (3) 

is the overlap matrix of the atomic basis, and the vectors 
0''(k)) are the dual of the atomic basis, that satisfy 

inm.m - (M^wm = s^^- (4) 

S measures the difference between the plane-wave eigen- 
states \ipa(k)) and their projection into the atomic-basis 
P(k)|V'Q(k)). More precisely, it gives the average of 
I (1 — P) IV'a(k)) I p over the eigenstates considered for the 
projection. The spilling S varies between and 1.5 = 
means that the LCAO projected wave-functions repro- 
duce the PW eigenfunctions exactly. 5 = 1 means that 
the basis is orthogonal to the Hamiltonian eigenstates. 

As mentioned above, the quality of an LCAO basis is 
evaluated not only for a particular system, but also for 
a particular set of Hamiltonian eigenstates. An LCAO 
basis will reproduce with very different accuracy the dif- 
ferent eigenstates of a Hamiltonian. The choice of eigen- 
states relevant to the definition of S depends on the par- 
ticular application. For ground-state properties of the 



system (total energy, geometry, electronic density, chem- 
ical analysis) the occupied eigenstates are the ones to be 
considered. For quasiparticle excitations or optical prop- 
erties, some of the lowest empty states should be also 
included in the projection. As shown later, the latter 
case typically requires a larger atomic basis. 

In addition to the analysis and possible selection of 
atomic basis sets, the methodology presented in this work 
allows for their systematic optimization in a straightfor- 
ward way. Given a type of atomic functions (say Slater- 
type orbitals, Gaussian-type orbitals) that depend on 
certain parameters (say exponents) the problem of op- 
timizing the basis becomes minimizing 5 as a function of 
those parameters. Note that the minimization requires 
evaluations of S for different values of the parameters, 
but for a single first-principles calculation. Taking into 
account that the evaluation of S represents a minimal 
computational effort compared with the actual PW cal- 
culation, it is obvious that this optimization procedure 
is much more convenient than the direct trial-and-error 
procedures using self-consitent LCAO codes. 

Using this procedure we have performed a systematic 
study comparing different kinds of atomic orbitals in dit. 
ferent solids. This will be presented in detail elsewhere.Ej 
In that analysis it is found that for many properties 
(spilling, charge density and populations, energy bands) 
the type of atomic orbitals that gives better results is the 
following: take the eigenfunctions of the free atoms us- 
ing the same basic apomximations as for the solid (LDA, 
same pseudopotentialO) and scale them down by scale 
factors A/ (one scale factor per different atomic orbital), 
that is, 0,(r) = Ni{Xi)q>f°"'{\ir), being 7V,(Az) the nor- 
malization factor. The optimum scale factors are deter- 
mined by minimization of S. Typically the orbitals com- 
press by a small amount except for the orbitals which are 
empty in the atomic calculation that compress substan- 
tially. These optimized orbitals are the ones used in the 
following, except stated otherwise. 

A further utility of the projection technique is its suit- 
ability for obtaining and analyzing LCAO Hamiltonians 
(tight-|binding parameters) and their associated energy 
bandsO The matrix elements of the projected Hamilto- 
nian 

H^^^^ik) = (cb,{k)\HP^\Mk)) (5) 

can be obtained using fast Fourier transform algorithms 
and other techniques of the PW method: 

H^^^'^ik) = E (0M(k)|k + G)(k + G|F^^|</..(k)) 

|k+G|2<£;„„^ 

(6) 

where G are reciprocal lattice vectors and Emax is an 
energy cut-off independent of the one used in the PW 
calculation, that must be large enough to guarantee a 
reliable representation of the atomic basis. It has to be 



2 



stressed that the proposed procedure obtains the Hamil- 
tonian matrix elements from first-principles, no free pa- 
rameters being fitted. 

The LCAO Hamiltonian matrix is obtained directly in 
Bloch space. By construction, it includes the information 
of the matrix elements between localized functions up to 
infinite neighbors. This is a very interesting feature of 
our method for the characterization of a basis in terms 
of its associated band structure: it conveniently separates 
the effect of the uncompleteness of the basis from the ef- 
fect of neglecting matrix elements beyond some scope of 
neighbors. These two sources of error are usually mixed 
together in the traditional ways of choosing atomic or- 
bitals from their LCAO bands. With our procedure each 
approximation can be analyzed separately. 

To obtain the real-space Hamiltonian matrix ele- 
ments (tight-binding parameters) from the reciprocal- 
space Hamiltonian, we perform an inverse Bloch trans- 
formation 

k 

where normalization factors which depend on the over- 
laps are omitted for clarity. The sum has to be extended 
to a sufficient number of k points, taking into account 
that the number of points depends on the real space range 
of the interactions. O For silicon and the ST0-4G basis 
used below, interactions are important up to third near- 
est neighbours and a few tens of k points have proven to 
be enough. 

In Figure 1 (a) and (b) we present the energy bands 
of silicon for an sp minimal basis and for an spd basis, 
respectively, with interactions up to infinite neighbors, 
as obtained from the projection. They are compared 
with the PW energy bands, that were calculated with 
a well converged plane- wave basis set. The atomic ba- 
sis was generated and optimized as explained above, the 
d orbital being obtained from an excited configuration 
of the silicon atom. The spilling of charge (considering 
the sum over occupied states in Eq. 1) is 5 = 0.0076 and 
S — 0.0007 for the sp and spd basis, respectively. The fig- 
ure shows that the sp basis reproduces the valence band 
much better than the lowest conduction bands, giving a 
rather poor description of the band gap. The inclusion 
of the d orbitals substantially improves both the valence 
band and the band gap, especially the latter. 

The atomic orbitals used above are relatively extended. 
For some purposes a reduction of their extension is de- 
sirable. Figure 2 shows the effects of a contraction of the 
basis. The energy of the X point rises, showing a ten- 
dency to establish a direct band gap, in addition to the 
expected narrowing of the bands and the overall increase 
of the gap. The shape and the ordering of the conduction 
bands is very sensitive to the contraction of the basis. See 
for example the change in ordering of the s and p con- 
duction states at F with the change in scale factor. The 
charge spillings are S = 0.0076 and S = 0.0292 for the 
optimum and the contracted basis, respectively. 



Even though the energy bands become worse with an 
overlocalization of the atomic orbitals, this can still be in- 
teresting if matrix elements for neighbors beyond a range 
are wanted to be negligible. The effect of this kind of 
approximation is shown in Figure 3, where the energy 
bands forj-silicon are presented for an ST0-4G sp mini- 
mal basisEEl with interactions up to first (a) , second (b) , 
and third (c) nearest neighbors. They result from the di- 
agonalization of the LCAO Hamiltonian matrix obtained 
in Eq. 5, and they are compared with the result of infinite 
neighbors. We have observed that the atomic basis ob- 
tained from the scaled eigenfunctions of the atom usually 
require a longer range of interactions. The use of ST0-4G 
orbitals reduces the range of non-negligible interactions. 
However, the quality of the bands (infinite neighbors) for 
the STO basis (dotted lines in Fig. 3) is worse than for 
the other basis, showed with solid lines in Fig. 1 (a), as 
compared with the plane-wave results, shown as dotted 
lines in Fig. 1. Notice that the optimized basis used for 
Figure 1 (a) correctly gives an indirect band gap, which 
is not usually obtained with standard minimal sp basis 
sets like the one used in Fig. 3. The charge spilling of 
the ST0-4G basis is S = 0.0513, considerably larger than 
for the basis used in Fig. 1 (a) {S = 0.0076). In LCAO 
calculations the extension of the atomic orbitals has to 
be chosen compromising between the quality required for 
the band-structure and valence properties, and the possi- 
bility of neglecting matrix elements beyond a range. The 
choice of basis depends on the particular application. 

Finally, we show the utility of the projection tech- 
nique for the analysis of the results of PW calculations 
by means of LCAO population analysis. We use the one 
proposed by Mulliken.la The analysis is performed on the 
occupied eigenstates projected into the subspace of the 
basis. It has to be done with care since the projected 
eigenstates |xa(k)) — P(k)\'ilja(k)) are not orthonormal 
(the deviation from orthonoriaaJity is related with the 
charge spilled in the projectiontS). Proper consideration 
of their associated overlap i?Q/3(k) = (xQ(k)|x/3(k)) is 
necessary to ensure charge conservation in the popula- 
tion analisys. This is accomplished by defining the den- 
sity operator as 

occ 

/5 = EEl:^"W)(x"(k)|, (8) 

k a 

where |x"(k)) = J2f3 ^paO^)\Xi3i'^)) represent the vec- 
tors of the dual set of the projected eigenstates. The 
LCAO density matrix is then written in terms of the 
dual LCAO basis: 

r,. = {r\pm, (9) 

the charge associated to an orbital /i then being 

Results of this procedure are shown in Table I, where 
charge transfers are calculated for different heteropolar 
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zincblende materials. They are compared with results of 
selfconsistent LCAO calculations showing a remarkable 
agreement. The basis functions used for these population 
analysis are the scaled atomic orbitals discussed above, 
which showed to be best suited for the purpose. 

In summary, it has been shown how the simple pro- 
jection of the eigenstates obtained from PW calculations 
into the space spanned by an LCAO basis can be use- 
ful for multiple purposes: (i) to evaluate and optimize 
LCAO basis sets, (ii) to analyze LCAO band structures, 
(iii) to obtain tight-binding parameters, and (iv) to per- 
form population analysis. The procedures are straight- 
forward and more systematic than the previous methods 
to perform equivalent tasks. The main idea of this paper 
could also be used for the improvement of first-principles 
algorithms. This possibility is presently being explored. 
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FIG. 1. Electronic band structure of silicon calculated 
projecting the plane- wave LDA Hamiltonian into an sp atomic 
basis (a), and into an spd atomic basis (b), both shown in solid 
lines. Dotted lines show the plane-wave LDA band structure. 
The basis functions are obtained as discussed in the text. 

FIG. 2. Effect of different contractions of an sp atomic 
basis in the silicon band structure. Solid lines are for opti- 
mum scale factors and dotted lines for a further contraction 
of the basis functions of a 1.2 factor. The basis functions are 
obtained as discussed in the text. 

FIG. 3. Band structure of silicon as a function of the 
range of non-neglected interactions: interactions up to first 
(a), second (b), and third (c) nearest neighbors. In dotted 
lines the result for infinite neighbors is shown as reference (in 
(c) they are hardly distinguishible from the solid lines). The 
atomic basis functions are standard STO-4G taken from Ref. 
10. 

TABLE I. Calculated charge transfer for some zincblende 
semiconductors. Qc and Qa stand for the valence charge on 
the cation and the anion, respectively, 5Q for the charge trans- 
fer with respect to neutral atoms, and S for the spilling. Num- 
bers in parenthesis were obtained from Hartree-Fock LCAO 
calculations,^^ except for GaAs, for which LDA LCAO was 
used.^^ 





Basis 


S 


Qc 


Qa 


SQ 


BN 


s,p 


0.0022 


2.19 (2.14) 


5.81 (5.86) 


0.81 (0.86) 


BP 


s,p 


0.0038 


3.51 (3.34) 


4.49 (4.66) 


-.51 (-.34) 


AlP 


s,p 


0.0035 


2.15 (2.20) 


5.85 (5.80) 


0.85 (0.80) 


SiC 


s,p 


0.0071 


2.30 (2.19) 


5.70 (5.81) 


1.70 (1.81) 


GaAs 


s,p 


0.0041 


2.58 


5.42 


0.42 




s,p,d 


0.0010 


2.78 (2.88) 


5.22 (5.12) 


0.22 (0.12) 
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